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Abstract. The outstanding capabilities of the Chandra 
X-ray observatory have greatly increased our potential to 
observe and analyze thermal radiation from the surfaces of 
neutron stars (NSs). Such observations allow one to mea- 
sure the surface temperatures and confront them with the 
predictions of the NS cooling models. Detection of gravita- 
tionally redshiftcd spectral lines can yield the NS mass-to- 
radius ratio. In rare cases when the distance is known, one 
can measure the NS radius, which is particularly impor- 
tant to constrain the equation of state of the superdense 
matter in the NS interiors. Finally, one can infer the chem- 
ical composition of the NS surface layers, which provides 
information about formation of NSs and their interaction 
with the environment. We will overview the recent Chan- 
dra results on the thermal radiation from various types 
of NSs -active pulsars, young radio-quiet neutron stars in 
supernova remnants, old radio-silent "dim" neutron stars- 
and discuss their implications. 



1. Introduction 

Observational study of thermal emission from isolated 
(non-accreting) neutron stars (NSs) became possible after 
the launch of Einstein (1979) and EXOSAT (1983), the 
first X-ray observatories equipped with X-ray telescopes. 
Contrary to nonthermal radiation of NSs, generated in the 
pulsar magnetospheres and observed from radio to 7-rays, 
thermal emission originates immediately at the surface, 
with the bulk of the energy flux in the soft X-ray band. 
Therefore, comparing the spectra and the light curves with 
the models for NS thermal radiation (see the contribution 
by Zavlin & Pavlov in these Proceedings; ZP02 hereafter), 
one can infer the NS surface temperature, magnetic field, 
gravitational acceleration, and chemical composition, as 
well as the NS mass and radius. The analysis of such ob- 
servations allows one to trace thermal evolution of NSs 
(Tsuruta 1998) and constrain the properties of the super- 
dense matter in the NS interiors (D.G. Yakovlev et al., 
these Proceedings). 

Observational manifestations of NSs are very diverse, 
and only some of these exotic objects are suitable for ob- 



serving the NS thermal emission. Most of the detected NSs 
are radio pulsars. Very young, active pulsars (with ages 
t < 10 3 yr — e.g., the Crab pulsar) are quite hot, with 
expected surface temperatures ~ 1-2 MK, but their non- 
thermal radiation is so bright that the thermal radiation 
is hardly observable. The nonthermal component in X-ray 
emission of middle-aged pulsars (r ~ 10 4 -10 6 yr — e.g., 
B0656+14, B1055-52) is much fainter, and their thermal 
radiation, with temperatures 0.3-1 MK, can dominate at 
soft X-ray and UV energies. The surfaces of old pulsars 
(t > 10 6 yr — e. g., B0950+08, B1929+10, J0437-4715) 
are too cold to be seen in X-rays, but their polar caps 
are expected to be heated up to a few million kelvins, hot 
enough to emit observable thermal X-ray radiation. 

In addition to active pulsars, a number of radio-quiet 
isolated NSs emitting thermal-like X-rays have been de- 
tected, with typical temperatures ~ 0.5-5 MK. They are 
usually subdivided in four classes: Anomalous X-ray Pul- 
sars (AXPs; see S. Mereghetti et al., these Proceedings), 
Soft Gamma-ray Repeaters (SGRs; Kouveliotou 1995), 
"dim" or "truly isolated" radio-silent NSs (i.e., not associ- 
ated with supernova remnants [SNRs]; Treves et al. 2000), 
and compact central sources (CCOs) in SNRs (Pavlov et 
al. 2002a) which have been identified with neither active 
pulsars nor AXPs/SGRs. Their observational manifesta- 
tions (particularly, multiwavelength spectra) are quite dif- 
ferent from those of "active" pulsars, and their properties 
have not been investigated as extensively, but the presense 
of the thermal component in their radiation provides a 
clue to understand the nature of these objects. 

First thermally emitting isolated NSs were detected 
with Einstein (e.g., Fahlman & Gregory 1981; Cheng & 
Helfand 1983; Helfand & Becker 1984; Seward, Charles, & 
Smale 1986; Cordova et al. 1989) and investigated with 
EXOSAT (Brinkmann & Ogelman 1987; Kellett et al. 
1987). In 1990's, ROSAT and ASCA detected X-rays from 
more than 30 radio pulsars, including a few thermal emit- 
ters (Ogelman 1995), and discovered several radio-silent 
NSs (see Becker & Pavlov 2002 for a review). Thermal 
emission from a few isolated NSs was also detected in the 
optical- UV range with the Hubble Space Telescope (e.g., 
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Pavlov, Stringfellow & Cordova 1996; Walter & Matthews 
1997). 

The new era in observing X-ray emission from NSs 
has started with the launch of two currently operating X- 
ray observatories of outstanding capabilities, Chandra and 
XMM-Newton. In this review we present recent results 
on thermal X-ray emission from isolated NSs of various 
types observed with Chandra. A more general overview of 
Chandra observations of NSs is presented in this volume 
by M.C. Weisskopf, while the XMM-Newton results on 
NSs are reviewed by W. Becker. 
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Fig. 2. Two-component (magnetic hydrogen NS atmo- 
sphere plus power law) model fit to the combined HRC- 
S /LETG and ACIS count rate spectra. The bottom panel 
shows the contributions from the thermal and nonthermal 
components (lilac and green lines, respectively). 



Fig. 1. Chandra HRC-I image of the Vela pulsar and 
surrounding nebula. 



2. Thermal emission from radio pulsars 

X-ray emission of rotation-powered pulsars generally con- 
sists of two components — thermal and nonthermal. The 
nonthermal component, generated by relativistic particles 
in the NS magnetosphere, is, as a rule, strongly pulsed 
and characterized by a power-law (PL) spectrum. Study- 
ing the magnetospheric component is important to un- 
derstand the mechanisms of the multiwavelength (optical 
through gamma-rays) pulsar radiation. The thermal com- 
ponent, emitted from the NS surface, can be observed in 
X-rays if it is not buried under the nonthermal compo- 
nent, as in very young pulsars, and if the temperature 
is not too low, as in old pulsars. Weak pulsations of the 
thermal emission can be caused by temperature nonuni- 



formities along the NS surface and/or anisotropy of the 
radiation intensities in a strong magnetic field (see ZP02). 
Stronger pulsations of thermal radiation can be observed 
in old pulsars which emit thermal X-ray radiation from 
small, hot polar caps around the magnetic poles. In this 
Section we present the results of our analysis of Chandra 
data on pulsars of different ages: the young Vela pulsar 
(t c = P/2P = 11 kyr), the middle-aged pulsars B0656+14 
(t c = HO kyr) and B1055-52 (r c = 500 kyr), and the very 
old millisecond pulsar J0437-4715 (r c = 4.9 Gyr). 

2.1. The Vela pulsar 

First Chandra observations of this famous 89-ms radio, 
optical and 7-ray pulsar in the Vela SNR (the distance 
d w 300 pc — Caraveo et al. 2001) were carried out with 
the imaging array of the High Resolution Camera (HRC- 
I) in January and February 2000. These observations re- 
vealed a very complicated structure of the PWN around 
the pulsar (see Fig. [l]) studied by Helfand, Halpern & Got- 
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Fig. 3. Multiwavelength energy spectrum of the Vela pul- 
sar. The solid line shows the NS atmosphere plus PL fit 
to the observed LETG and ACIS spectra. The dotted line 
corresponds to same model spectrum corrected for the in- 
terstellar absorption and extrapolated to lower and higher 
energies. The dash-dot lines show the extrapolated opti- 
cal and EUV absorbed spectra (see Pavlov et al. 2001b for 
details and references). 
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Fig. 4. Light curve of the Vela pulsar from two HRC- 
I observations of January and February 2000. The zero 
phase corresponds to the radio peak. 

thelf (2001) and Pavlov et al. (2000, 2001a). The superb 
spatial resolution of Chandra made it possible, for the first 
time, to resolve the pulsar from its surroundings and in- 
vestigate the pulsar's radiation. 

The spectral data on the Vela pulsar were obtained in 
two observations (Pavlov et al. 2001b). A 25 ks HRC-S ob- 
servation with Low Energy Transmission Grating (LETG) 
of January 2000 provided the high-resolution pulsar spec- 
trum in the energy range E = 0.25-2.0 keV. The anal- 
ysis of the LETG spectrum (« 7, 000 source counts) did 



not reveal statistically significant spectral lines. The spec- 
trum is best fitted with thermal models (blackbody, hy- 
drogen NS atmosphere), but the observed spectrum some- 
what exceeds the model spectra at E > 1.5 keV (see upper 
panel in Fig. |^), indicating the presence of a second com- 
ponent with a harder spectrum. The spectrum at higher 
photon energies was provided by a 37 ks observation (Oc- 
tober 1999) with the Advanced CCD Imaging Spectrome- 
ter (ACIS), taken in the Continuous Clocking (CC) mode 
that allows timing at the expense of imaging capability in 
one dimension. The 1-D count distributions of the pulsar 
and the nebula show that the pulsar is much brighter than 
the PWN background at E < 1.5 keV, and it is discernible 
up to 8 keV, where the charged particle background takes 
over (see Fig. 1 in Pavlov et al. 2001b). The pulsar spec- 
trum in the 0.5-8 keV band (« 9, 000 source counts) does 
not fit one-component models, but it fits fairly well two- 
component models (e.g., a PL plus a thermal component 
— see the middle panel in Fig. 2). 

The analysis of the combined (LETG plus ACIS) 
spectra proves that the pulsar's soft X-ray emission, at 
E < 1.8 keV, is indeed dominated by a thermal compo- 
nent (see bottom panel in Fig. ||). The parameters of 
this component are drastically different for the blackbody 
fit (T b °g = 1.4-1.5 MK, i?g£ = [2.3-2.9] d 3Q0 km) and 
the NS hydrogen atmosphere fit (TJ? = 0.65-0.71 MK, 
R°° — [17-20]d3oo km). Such strong difference is explained 
by the fact that the spectrum of a hydrogen atmosphere 
decreases with increasing photon energy slower than the 
ideal Wien spectrum because the radiation at higher ener- 
gies comes from hotter layers. As a result, the blackbody 
fit to a hydrogen atmosphere spectrum yields a temper- 
ature higher than the true effective temperature and an 
area smaller than the true emitting area (see ZP02 for 
more details). 

The small blackbody radius might hint that the ob- 
served radiation is emitted from hot spots (e.g., polar 
caps heated by relativistic particles produced in the pulsar 
magnetosphere) , whereas the radius obtained in the hy- 
drogen atmosphere fit implies radiation emitted from the 
entire NS surface. Rigorously speaking, the surface (at- 
mosphere) of a NS is not a black body, but a blackbody 
spectrum could mimic the spectrum of a heavy-element 
(e.g., iron) atmosphere if it is observed with low spectral 
resolution. However, the lack of significant spectral lines in 
the LETG spectrum indicates that there are no heavy ele- 
ments, neither He nor metals, in the radiating atmospheric 
layers. On the other hand, even strongly magnetized hy- 
drogen does not have spectral features in the investigated 
energy band (all the hydrogen features are below 0.17 keV 
for the field of B = 3 x 10 12 G, adopted for the Vela pul- 
sar). Therefore, the hydrogen atmosphere interpretation 
looks more plausible. With this interpretation, the effec- 
tive temperature of the Vela pulsar is well below the pre- 
dictions of the "standard" models of the NS cooling (e.g., 
Tsuruta 1998). 
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An important result of this analysis is the detec- 
tion of the hard spectral component which dominates at 
E > 1.8 keV. The slope (photon index) of this PL compo- 
nent is 7 = 1.2-1.8 if the NS hydrogen atmosphere model 
is adopted for the soft component. With this interpreta- 
tion, the overall X-ray spectrum of the the Vela pulsar 
is quite similar to those of the three middle-aged pulsars, 
Geminga, B1055-52 and B0656+14 (see below), where the 
thermal component was detected, although the ratio of the 
nonthermal X-ray luminosity to the spin-down energy loss 
rate, L x /E ~ 3 x 10 -6 in the 0.2-8.0 keV range, is a fac- 
tor of 10 3 smaller than observed for typical X-ray detected 
pulsars (Becker & Triimper 1997). If the L x /E ratio were 
as high as for other pulsars, the thermal radiation would 
be undetectable in the phase-integrated spectrum. 

It is illuminating to compare the nonthermal X-ray 
spectrum of the pulsar with those in other energy ranges. 
Figure [| shows the spectral flux of the pulsar from op- 
tical to gamma-rays. Remarkably, the extrapolation of 
the PL spectrum obtained in the atmosphere-plus-PL fit 
matches fairly well with the optical and hard- X-ray (soft- 
7-ray) time-averaged fluxes. Moreover, under the hypothe- 
sis that the PL spectrum has the same slope in the optical 
through X-ray range, we can constrain the photon index 
very tightly: 7 = 1.35-1.45. 

The multicomponent structure of the pulsar radiation 
is also seen from its pulse profile. The HRC-I light curve is 
shown in Figure [| (see also Helfand et al. 2001; Sanwal et 
al. 2002a). It has at least three peaks — quite unusual for 
a pulsar, but not so surprising for the Vela pulsar which 
shows up to four peaks in the hard X-ray bands observed 
with RXTE (Harding et al. 2002). 

So far, only the analysis of phase-integrated spec- 
tra and energy-integrated light curves has been possible. 
Phase-resolved spectral analysis and energy-resolved tim- 
ing would be of crucial importance to elucidate the com- 
plicated emission mechanism(s) of both the nonthermal 
and thermal components. In particular, the analysis of the 
phase-energy dependence of the thermal emission, based 
on the NS atmosphere models can make it possible to de- 
termine the surface distributions of the temperature and 
the magnetic field. 

2.2. Middle-aged pulsars B0656+U and B1055-52 

Among about 40 rotation-powered pulsars detected with 
the previous and currently operating X-ray missions, 
middle-aged pulsars (r ~ 0.1-1 Myr) are of special in- 
terest because their thermal emission is expected to dom- 
inate in soft X-rays. The best studied objects of this class 
are the "Three Musketeers" Q: PSR B0656+14, Geminga, 
and PSR B1055-52. 



1 As dubbed by J. Triimper 
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Fig. 5. Three-component, TS+TH+PL, fit to the com- 
bined LETG and ACIS-S spectra of PSR B0656+14. The 
red curves show the model countrate spectra for the two 
instruments. 

2.2.1. PSR B0656+14 

PSR B0656+14 (r c = 110 kyr) is the brightest of the 
Three Musketeers. Observations with Einstein, ROSAT, 
and ASCA have shown that its X-ray spectrum consists of 
at least two components. The thermal component (E < 1 
keV) resembles a soft blackbody spectrum with an appar- 
ent temperature of about 0.9 MK. Because of the narrow 
energy band of ROSAT and poor angular resolution of 
ASCA, the shape of the hard component (at E > 1 keV) 
could hardly be determined from those data. For instance, 
the ASCA spectra above 1 keV can be fitted equally well 
with a hard blackbody with temperature of about 2.5-4.0 
MK or a PL with photon index 7 = 1.3-3.0. IR-optical-UV 
observations of the pulsar (Pavlov et al. 1997; Koptsevich 
et al. 2001) revealed a nonthermal PL component with 
7 = 1.4-1.5, presumably emitted from the pulsar magne- 
tosphere. 

To search for spectral features in the X-ray spectrum 
of PSR B0656+14 and investigate the shape of its high- 
energy tail, the pulsar was observed with HRC-S/LETG 
in November 1999 (38 ks exposure) and ACIS-S (in CC 
mode) in December 2001 (25 ks exposure). The analy- 
sis of these data showed no significant features in the 
source spectrum. To obtain a good fit to the combined X- 
ray spectrum in the 0.2-6 keV band, a three-component 
model is required. The high-energy tail (E > 2 keV) fits 
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Fig. 6. Spectral flux for the TS+TH+PL model for PSR 
B0656+14 extrapolated to the optical range. Blue and red 
curves show the absorbed and unabsorbed spectra, respec- 
tively. Crosses are the IR-optical-UV fluxes (Koptsevich et 
al. 2001). 



best with a PL spectrum with a slope close to that of 
the optical spectrum. To constrain it more tightly, we 
assumed that the slope does not change from optical 
to X-rays, included the optical data in the fit, and ob- 
tained 7 = 1.5. The low-energy part of the spectrum 
(E < 0.7 keV) fits well with a blackbody spectrum with 
temperature T s °° = 0.85 MK and equivalent radius of 
emitting sphere R™ = 12.0 (d/500 pc) km. This ther- 
mal soft (TS) component is presumably emitted from the 
whole visible NS surface. To fit the spectrum at inter- 
mediate energies, one more component is required (see 
Fig. ||). Very likely, this is a thermal hard (TH) com- 
ponent with temperature T£° ~ 1.65 MK (as given by 
the blackbody model) emitted from pulsar's polar caps 
of radius R pc = ~ 1.0 (<i/500 pc) km heated by rela- 
tivistic particles streaming down from the pulsar acceler- 
ation zones. The bolometric luminosity of the TS com- 
ponent is Lboi.s — 5.5 x 10 32 (<i/500pc) 2 erg s" 1 , vs. 
L bo i )h ^ 4.9 x 10 31 (d/500pc) 2 erg s" 1 for the TH com- 
ponent. The luminosity of the nonthermal component is 
L x ~ 5.3 x 10 30 (d/500pc) 2 erg s" 1 in the 0.1-10.0 keV 
range. This three-component model fits the data from IR 
through X-ray energies (Fig. ^|). 

Using the magnetic hydrogen atmosphere models 
(ZP02), instead of blackbodies, for the thermal compo- 
nents gives formally acceptable fits, with effective tem- 
peratures lower by a factor of about 2 (e.g., ~ 0.4-0.5 
MK for the TS component). However, these fits imply ei- 
ther very small distances, ~ 70-130 pc for R = 10 km, or 
very large radii, R s ~ 40-70 km for the TS component, 
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Fig. 7. X-ray pulse profiles for PSR B0656+14 in different 
energy bands (four upper panels). The radio pulse pro- 
file at 1.642 GHz (Gould & Lyne 1998) demonstrates the 
phase difference between the X-ray and radio peaks. The 
X-ray data are from the Chandra ACIS observation in the 
CC mode (December 2001). 



at d = 500 pc. Therefore, applicability of these models to 
PSR 0656+14 is questionable. 

The multicomponent structure of the pulsar X-ray ra- 
diation is also seen from the behavior of the pulse profiles 
shown in Figure |?]. The shape of the light curves is energy 
dependent, with a significant phase shift (~ 0.2 of the pe- 
riod) at 0.5-0.7 keV. The source pulsed fraction f p shows 
a complex behavior: it decreases with increasing energy 
from about 20% at E ~ 0.1-0.3 keV down to about 7% 
at E around 1 keV and increases again to / p ~ 35% at 
E > 1 keV. The X-ray peak lags the radio peak by 0.2-0.4 
of the period, depending on energy (Fig. 0). 

To establish the origin of the peak and the parame- 
ters of the X-ray radiation, one should analyze the spec- 
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Fig. 8. Image of a 10" x 10" area around PSR B0656+14 
from the 5 ks ACIS observation in the Timed Exposure 
mode (December 2001). 
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Fig. 9. Three-component model fit to the Chandra ACIS 
and ROSAT PSPC spectra of PSR B1055-52. 



tra and the light curves together, with the aid of realistic 
models for X-ray emission. However, the spectra of the 
pulsar observed with the pre- Chandra X-ray observato- 
ries, and even the Chandra ACIS-S spectrum observed in 




log E (keV) 

Fig. 10. Spectral flux from PSR 1055-52 for the 
TS+TH+PL model extrapolated to the optical and 7-ray 
ranges, together with the data from HST, OSSE, COMP- 
TEL and EGRET. The blue and the red curves give the 
absorbed and unabsorbed spectra, respectively. 



CC mode, potentially have a contaminating component — 
synchrotron radiation from a PWN of a small angular size 
around the pulsar, which could not be resolved in those 
observations. 

A hint on extended emission around PSR 0656+14 is 
seen in the zero-order LETG/HRC-S image (Marshall & 
Schulz 2002), but that image is strongly contaminated by 
the diffraction spikes. The pulsar was also observed in a 
5 ks exposure with ACIS-S3 in Timed Exposure mode 
(December 2001) providing the data for spatial analy- 
sis of the X-ray emission in the vicinity of the pulsar. 
The analysis reveals a 3"-4" extended structure around 
the piled-up image of the pulsar — first strong indica- 
tion on PWN around a middle-aged pulsar (see Fig. ||). 
The pulsar's image, confined within a small, slightly elon- 
gated area of about 1" diameter (at the center of Fig. 
||) is surrounded by a putative PWN of rather compli- 
cated morphology. The estimated luminosity of the PWN 
is L X:P wn - 2 x 10 31 (d/500 pc) 2 erg s" 1 in the 0.1-10.0 keV 
range. A considerable fraction of the PWN radiation can 
contaminate the pulsar's radiation observed with the pre- 
Chandra observatories and even with ACIS in CC mode. 
Therefore, to unambiguously separate the PWN contribu- 
tion, one should observe the pulsar and its surroundings in 
high-resolution imaging mode for longer time. The above 
estimates of the pulsar spectral parameters should be con- 
sidered with caution until the properties of the putative 
PWN are studied and its radiation is separated from that 
of the pulsar. 



G.G. Pavlov, V.E. Zavlin & D. Sanwal: Thermal Radiation from Neutron Stars: Chandra Results 



7 




0.5 1 1.5 0.5 1 1.5 

Radio Phase Radio Phase 



Fig. 11. Energy-resolved light curves of PSR B1055-52 
plotted vs. radio phase. 

2.2.2. PSR B1055-52 

Analysis of the Chandra data on PSR B1055-52, taken 
in January 2000 with ACIS-S in CC mode (42 ks expo- 
sure), yields results very similar to those obtained for PSR 
B0656+14 (see Sanwal et al. 2002c for details). First, the 
ACIS spectrum, combined with that obtained with the 
ROSAT PSPC, requires a three-component model (Fig. 
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Fig. 12. ACIS image of PSR 0437-4715. The bright strip 
is the trailed image of a nearby AGN. The dashed curve 
marks the leading edge of the H a bow-shock. The arrow 
shows the direction of the pulsar's proper motion. 



^) . The high-energy tail of the spectrum (E > 2 keV) is 
best described with a PL with 7 = 1.66; such a compo- 
nent smoothly connects optical (HST FOC) and gamma- 
ray (CGRO EGRET) data (Fig. |l^). The spectrum at 
lower energies fits well with two thermal (blackbody) com- 
ponents, soft and hard. The best-fit model parameters to 
the phase-integrated spectrum are T s °° = 0.84 MK and 
= 12.0 (d/lkpc) km for the soft thermal component, 
T h °° = 1.62 MK and Pg° = 0.8 (d/1 kpc) km for the hard 
thermal component, hh = 1.3 x 10 20 cm~ 2 . 

The energy-resolved light curves of the pulsar derived 
from the ACIS data reveal energy-dependent variations in 
the pulsed fraction, the phase of the pulse, and the pulse 
shape (Fig. [ll]) . In particular, the pulsed fraction increases 
from - 15% at about 0.3-0.4 keV to - 50% at > 1.0 keV. 
For this pulsar, the X-ray pulse at E > 0.3 keV is in phase 
with the main radio pulse, contrary to PSR B0656+14. 

The phase-resolved spectroscopy shows small varia- 
tions in the parameters of the soft thermal component, 
whereas the emitting area for the hard thermal component 
varies significantly with the rotational phase changes 
between about 0.5 and 1.5 km), which can be interpreted 
as changing projection of the hot spot area. 

Thus, we see remarkable similarity in the X-ray proper- 
ties of the two middle-aged pulsars, B0656+14 and B1055- 
52. Interestingly, the temperature of the soft thermal com- 
ponent is approximately the same for these pulsars, al- 
though the characteristic age, r c = P/2P, is a factor of 
5 larger for B1055-52. A possible explanation of this fact 
is that these NSs have different properties (e.g., the mass 
of B1055-52 is lower — see D.G. Yakovlev et al., this vol- 
ume). An alternative explanation is that the true age of 
a pulsar can be quite different from its characteristic age 
(see below), so that one should be very cautious comparing 
the NS cooling models with the observations of pulsars. 

2.3. Millisecond pulsar 30^31-^115 

Millisecond (recycled) radio pulsars are distinguished from 
ordinary pulsars by their very short and stable periods, 
P ^ 10 ms, P ~ 10~ 20 s s _1 . They are presumably very 
old objects, with spin-down ages r ~ 10 9 -10 10 yr and 
low magnetic fields, B ~ 10 8 -10 10 G. Similar to ordinary 
pulsars, a millisecond pulsar can emit nonthermal X-rays 
from its magnetosphere. Additionally, thermal X-rays can 
be emitted from pulsars' polar caps heated up to X-ray 
temperatures by relativistic particles generated in the pul- 
sar acceleration zones. 

At a distance d ~ 140 pc, PSR J0437-4715 (P = 5.75 
ms) is the nearest and brightest millisecond pulsar known 
at both radio and X-ray wavelengths. It is in a 5.74 d 
binary orbit with a low-mass white dwarf companion. Op- 
tical observations in H a have revealed a bow-shock PWN, 
caused by the supersonic motion of the pulsar through the 
interstellar medium, with the bow-shock apex at about 
10" south-east of the pulsar, in the direction of the pul- 
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Fig. 13. Two-temperature hydrogen polar cap model fit 
to the combined ACIS and PSPC data on PSR J0437- 
4715. 
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Fig. 14. Power-law plus two-temperature hydrogen polar 
cap model fit to the combined ACIS and PSPC spectrum 
of PSR J0437-4715. 



sar's proper motion. Its X-ray radiation was discovered 
with ROSAT (Becker & Triimper 1993). The pulsar was 
also observed with EUVE (Halpern, Martin, & Marshall 
1996). The pre- Chandra data, obtained in a narrow energy 
range of 0.1-2.0 keV, can be well fitted with a PL with 
photon index 7 = 2.2-2.5, or with a polar cap model with 
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Fig. 15. Pulse profiles of PSR 0437-4715 from the Chan- 
dra and radio observations. 



a nonuniform temperature distribution (Zavlin & Pavlov 
1998 [ZP98]). To discriminate between these two compet- 
ing models, the pulsar was observed with Chandra ACIS 
(May 2000; 26 ks exposure; see Zavlin et al. 2002). 

The source spectrum prevails over background at en- 
ergies up to 7 keV and reveals no significant spectral fea- 
tures. No structures which could be associated with the 
bow shock observed in Hq, are found (see Fig. 12). 



The combined analysis of the ACIS and ROSAT PSPC 
data rules out a single PL model (min xt — 3.8), indicat- 
ing that more components are needed to get a satisfactory 
combined fit. 

A two-component model, PL plus two polar caps of 
uniform temperature, provides an acceptable fit to the 
combined data. However, comparison with the EUVE data 
shows that this model is only marginally acceptable. In 
addition, the inferred values of hydrogen column density 
exceed the values nn = (1 _ 3) x 10 19 cm -2 derived from 
independent measurements of the interstellar absorption 
towards objects in the vicinity of the pulsar (see ZP98 for 
more details). 

ZP98 suggested that the temperature distribution in 
a polar cap should be nonuniform because the heat, re- 
leased by decelerating magnetospheric particles in sub- 
photospheric layers, propagates along the surface out from 
the small polar areas, where the energy is deposited, and 
heats a larger area. They used a two-zone approximation 
for the polar cap temperature distribution, a hot "core" 
plus a colder "rim" . When applied to the ACIS data, this 
model shows that the observed spectrum significantly ex- 
ceeds the model at E > 2 keV, i.e., one more component 
is required (Fig. |l3|). A natural choice for the additional 
component is a PL. Adding the nonthermal component 
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Fig. 16. Light curves of CXO J1210-5226 from two Chan- 
dra observations. 



greatly improves the fit (Fig. [lj]). At fixed n H = 2x 10 19 
cm~ 2 , the best fitting parameters are T coro = 2.1 MK, 
T rim = 0.54 MK, R COIC = 0.12 km, R lim = 2.0 km, and 
7 = 2.2. The bolomctric luminosity of two polar caps is 
Lboi = 2.3 x 10 30 erg s" 1 - 0.6 x 10" 3 £;, and the PL lumi- 
nosity in the 0.1-10 keV band is L x = 0.7 x 10 30 erg s -1 
~ 0.2 x lO -3 ^. This model is in good agreement with the 
EUVE data. 

In the polar cap model with temperature decreasing 
outwards from the cap center, the thermal component 
dominates between 0.1 keV and 2.5 keV, providing about 
75% of the X-ray flux, while the PL component dominates 
outside this band. An important point here is that this 
model requires the thermal radiation to be emitted from 
a hydrogen (or helium) NS atmosphere. The high spec- 
tral resolution of the ACIS data rules out an atmosphere 
comprised of heavier chemical elements. 

The pulsar was also observed with the HRC-S detector 
(October 2000; 20 ks exposure) in the timing mode that 
provides accurate timing results. Figure |l5| presents the 
HRC-S and radio pulse profiles. The phases of the X-ray 
and radio pulses are found to be virtually the same — their 
difference (0.3% of the pulsar's period) is smaller than the 
absolute timing accuracy of Chandra. This implies that, 
if the main contribution to the HRC-S band is due to 
the thermal polar cap radiation, then the radio emission 
must be generated close to the NS surface — e.g., the time 
difference of < 0.1 ms between the X-ray and radio phases 
corresponds to a distance of < 30 km, much smaller than 
the light cylinder radius of 275 km. 



3. Compact central objects in supernova 
remnants 

Chandra has observed five CCOs. Their most important 
characteristics are listed in Table 1. Although some prop- 
erties of CCOs (in particular, their spectra) are very sim- 
ilar to each other, it does not necessary mean that they 
represent a uniform class of objects. The most outstand- 
ing is J1617-5102 in RCW 103, with its highly variable 
flux (time scales hours to months), putative 6.6 hr period, 
and a possible optical counterpart (J ~ 22.3). One can 
presume that this source is not a truly isolated NS (i.e., 
it is powered by accretion, at least in its high state), al- 
though it is not a standard accretion-powered pulsar. The 
other four CCOs have shown neither long-term variabili- 
ties nor indications of binarity. If, based on the observed 
spectra, we adopt a plausible hypothesis that their radi- 
ation is thermal and assume that it can be adequately 
described by the same spectral model (e.g., BB or hydro- 
gen NS atmosphere), then we have to conclude that the 
emitting area is growing with age. If the radiation emerges 
from a hydrogen atmosphere, then the effective radius is 
consistent with a NS radius for two oldest CCOs, J0821- 
4300 and J1210-5226, being smaller for two youngest ones, 
J2323+5848 and J0852-4617. One might assume that the 
latter two are black holes, not NSs, but then one would 
have to invoke accretion as the energy source, which does 
not look consistent with the lack of variability. 

A recent review on CCOs, including the Chandra re- 
sults, has been presented by Pavlov et al. (2002a). There- 
fore, here we will describe only the most recent results on 
the best-investigated CCO, J1210-5226 (a.k.a. IE 1207.4- 
5209), the first isolated NS which shows spectral lines. 

3.1. J 1210-5226: the puzzling pulsar m G 296. 5- 10 

This radio-quiet X-ray source at the center of G296. 5+10.0 
(a.k.a. PKS 1209-51/52) was discovered by Helfand & 
Becker (1984). Mereghetti, Bignami, & Caraveo (1996) 
and Vasisht et al. (1997) fitted the ROSAT and ASCA 
spectra of J1210-5226 with a blackbody model of T{£ ~ 
3.0 MK from an area with radius R^ ~ 1.6 (d/2 kpc) km. 
Zavlin, Pavlov & Triimper (1998) interpreted the ob- 
served spectra as emitted from a light-element (hydro- 
gen or helium) atmosphere. For a NS of mass 1.4 Mq and 
radius 10 km, they obtained a NS surface temperature 
T°£ = (1.4-1.9) MK and a distance of 1.6-3.3 kpc, consis- 
tent with the distance to the SNR, 2.1+J' g kpc (Giacani 
et al. 2000). 

Zavlin et al. (2000) observed J1210-5226 with Chan- 
dra ACIS-S in CC mode in January 2000 (32 ks exposure) 
and discovered a period of about 424 ms, which proved 
that the source is a NS. Second Chandra observation with 
the same observational setup (January 2002; 30 ks) pro- 
vided an estimate of the period derivative, P ~ (0.7- 
3) x 10" 14 s s" 1 (Pavlov et al. 2002b). This estimate im- 
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Object 


Host SNR 


Age 


Period 




poo 


-^bol.bb 


J cS 




r oo 
^bol 






kyr 




MK 


km 


10 33 erg s _1 


MK 


km 


10 33 erg s" 1 


J2323+5848 


Cas A 


0.32 




5.7 


0.5 


1.6 


3.5 


1 


2 


J0852-4617 


G266.1-1.2 


~1 


301 ms [?] 


4.6 


0.3 


0.3 


3.1 


1.5 


0.8 


J1617-5102 


RCW 103 


1-3 


6.6 hr [?] 


4.6-7.0 


0.2-1.6 


0.5-30 


3.5 


1-8 


1-60 


J082 1-4300 


Pup A 


1-3 




4.4 


1.4 


4.2 


2.0 


10 


8 


J1210-5226 


G296.5+10.0 


3-20 


424 ms 


2.9 


1.6 


1.3 


1.6 


11 


1 



Table 1. Compact central objects in SNRs with best-fit parameters for blackbody (columns 5-7) and hydrogen NS 
atmosphere models (columns 8-10). The periods marked with "[?]" require confirmation. More details on these objects 
can be found in Pavlov et al. (2002a). 
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Fig. 17. Left: Fit to the Chandra spectrum of J1210-5226 with a blackbody model. Middle: The same model with 
two absorption lines added (at energies of 0.7 keV and 1.4 keV). Notice the significant improvement in the fit quality. 
Right: The two- line model for the photon flux (see text). 



plies that the characteristic age of the NS, r c ~ 200-1600 
kyr, is much larger than the 3-20 kyr age of the SNR 
(Roger et al. 1998), while the "canonical" magnetic field, 
B = 3.2 x 10 W (PP) 1 / 2 G = (2-4) x 10 12 G, is typical for 
a radio pulsar. These observations revealed smooth pulsa- 
tions of the X-ray flux from J1210-5226, with pulse shape 
depending on photon energy (Fig. |l6| ). 

The analysis of the ACIS spectra (Sanwal et al. 2002b) 
shows that continuum models (e.g., blackbody, power-law, 
fully ionized NS atmosphere) fail to fit the source spec- 
trum, which deviates very significantly from any of the 
models, revealing two absorption features near 0.7 keV 
and 1.4 keV (Fig. [l7]). Including absorption lines in the 
spectral model improves the quality of the fit substan- 
tially. Assuming Gaussian profiles for the absorption co- 
efficients, we obtain the central energies 0.73 ± 0.01 and 
1.40 ± 0.01 keV, equivalent widths 108 ± 15 and 112 ± 24 
eV, and central optical depths 0.40 ±0.06 and 0.43 ±0.06, 
respectively. 



Attempts to explain the absorption features as caused 
by the intervening interstellar (or circumstellar) material 
lead to huge overabundance for some elements. Therefore, 
the observed lines are most likely intrinsic to the NS. 

One could assume that the observed lines are cyclotron 
lines produced in a strongly ionized NS atmosphere. If one 
assumes these are electron cyclotron lines, the features 
might be interpreted as the fundamental and the first har- 
monic of the electron cyclotron energy E ce — 1.16-Bn keV 
in a magnetic field B n = S/(10 n G) - 0.6 (1 + z) (z is the 
neutron star rcdshift parameter), or as two fundamentals 
emitted from two regions with different magnetic fields, 
Bn = 0.6 (1 + z) and 1.2 (1 + z), broadened by the radia- 
tive transfer effects and/or nonuniformity of the magnetic 
field. However, it is difficult to reconcile this interpretation 
with the expected strength of the surface magnetic field 
B > 3 x 10 12 G. In addition, the oscillator strength of the 
first harmonic is smaller than that of the fundamental by 
a factor of ~ E ce /(m e c 2 ) ~ 2 x 10~ 3 (at E ce > kT), so 
that it is hard to explain why the 1.4 keV feature is as 



G.G. Pavlov, V.E. Zavlin & D. Sanwal: Thermal Radiation from Neutron Stars: Chandra Results 



11 



Object 


J bb 


h a 




Period 


B 


Ref. 




MK 




1n 2l] — 2 

10 cm 


s 


mag 




RX J 1856. 5-3754 


0.7 


14 


1.3 




25.8 


1 


RX J0720.4-3125 


0.9 


11 


1.1 


8.37 


26.6 


2,3 


RX J1605.3+3249 


1.1 


9 


1.1 




> 25 


4 


RX J0806.4-4123 


0.9 


3 


2.5 


11.37[?] 


> 24 


5,6 


RX J0420.0-5022 


0.7 


1 


1.7 


22.69[?j 


> 25.2 


7 


RX J1308.8+2127 


1.4 


5 


2.1 


5.16 


> 26 


8,9 


RX J2143.0+0654 


1.1 


9 


4.6 




> 22. 8 b 


10 



Table 2. Main properties of the seven "dim" NSs discovered with ROSAT. The spectral parameters are obtained 

from blackbody fits to the ROSAT PSPC data. The periods marked with "[?]" require confirmation. 

Refs.: 1 - Walter et al. (1996); 2 - Haberl et al. (1997); 3 - Kulkarni & van Kerkwijk (1998); 4 - Motch et al. (1999); 5 

- Haberl et al. (1998); 6 - Haberl & Zavlin (2002); 7 - Haberl et al. (1999); 8 - Schwope et al. (1999); 9 - Hambaryan 

et al. (2001); 10 - Zampicri et al. (2001). 

a Absorbed flux in 0.1-2.0 keV, in units of 10~ 12 erg cm~ 2 s _1 . 

b R magnitude. 



strong as the 0.7 keV feature if we assume the two lines 
are associated with the same magnetic field. 

One can assume that the spectral features are asso- 
ciated with ion cyclotron energies, E C i = 0.63(Z/yl)£>i4 
keV, where Z and A are the ion's charge and mass num- 
bers. The surface magnetic field needed for this interpre- 
tation is <; 10 14 G, much larger than the conventional 
magnetic field inferred for a centered dipole from the P, 
P measurements. Such a strong surface magnetic field can 
be consistent with the magnetic moment /i = 3 x 10 30 G 
cm 3 if the magnetic dipole is off-centered to about 3 km 
below the surface. However, the ion-cyclotron interpreta- 
tion of the observed features is not straightforward even 
in this case (see Sanwal et al. 2002b for more details). 

Thus, the most viable possibility is that the observed 
lines arc atomic lines formed in the NS atmosphere. The 
features cannot be explained as emerging from a hydrogen 
atmosphere because, at any magnetic field and any rea- 
sonable gravitational redshift, there is not a pair of strong 
hydrogen spectral lines whose energies would match the 
observed ones. Pavlov et al. (2002c) suggest that these are 
the strongest lines of once-ionized helium in a superstrong 
magnetic field, (1.4-1.7) x 10 14 G. In this interpretation, 
the observed line energies correspond to the gravitational 
redshift z = 0.12-0.23, which corresponds to R/M = 8.8- 
14.2 km Mq 1 . To confirm this result, deep observations 
with high spectral resolution are needed. 

4. "Dim" isolated neutron stars 

ROSAT discovered seven objects which are often called 
"dim" or "truly isolated" NSs. These sources, listed in 
Table 2, are characterized by blackbody- like spectra with 
temperatures around 1 MK. Low values of hydrogen col- 
umn densities, nn ~ 10 20 cm~ 2 , suggest small distances 
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io- 



to these objects. Their X-ray fluxes, / x ~ 10" 
erg s"" 1 cm" 2 , correspond to rather low (bolometric) lu- 
minosities, L <~ 10 30 — 10 31 erg s _1 for a fiducial distance of 
100 pc. The fluxes do not show variations on either short 
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Fig. 18. Chandra LETG/HRC-S spectrum of 
RX J1856. 5-3754 collected in a 500 ks exposure. 
The red curves represent the best blackbody fit to the 
observational data (the model parameters are given in 
the bottom panel). 



(minutes-hours) or long (months-years) time scales, and 
they are much higher than the optical fluxes, / x // op t >=S 
10 3 . Such properties strongly suggest that these objects 
are isolated, nearby NSs. Rotational periods in J0720.4- 
3125 and J1308.8+2127 were discovered with ROSAT and 
Chandra, respectively (see refs. in Table 2). There are also 
indications on periodicities in J0420. 0-5022 and J0806.4- 
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4123, although those periods are to be confirmed yet. The 
source powering the thermal radiation of these NSs has 
been a matter of debate (e.g., Treves et al. 2000). Two 
main options under consideration are the internal heat of 
a cooling NS and accretion from the interstellar medium 
(ISM). A standard estimate on the luminosity powered 
by the Bondi-Hoyle accretion, L x ~2x 10 31 ni; u0 3 erg s _1 
(here n is the number density of the ISM in cm -3 , viq is 
the relative NS velocity in units of 10 km s _1 ), requires 
rather low NS velocities. Recent estimates on the distance 
to J1856. 4-3754 and its proper motion (see below) imply 
a transverse velocity v± ~ 185 km s _1 . Therefore, accre- 
tion cannot be the energy source for this object unless 
the star moves in a surprisingly dense medium (a cloud?) 



with 



10 cm . First estimates of period derivatives 



for J0720.4-3125 (Zane et al. 2002; Kaplan et al. 2002b) 
and J1308.8+2127 (Hambaryan et al. 2002) imply surface 
magnetic fields^ ~ 10 13 — 10 14 G, which favors models in- 
voking young (or middle-aged) cooling NSs rather than 
accreting NSs. 

Below we present results from Chandra observations 
of the most famous (and best-investigated) objects of this 
class, J1856.5-3754 and J0720.4-3125. 

4.1. RX J '1856.5-3754 

This object (we will call it J1856 throughout this Sec- 
tion) was first observed with Chandra in March 2000 in 



We note that these results are preliminary and require 
further confirmation. 
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Fig. 19. Fits to the LETG spectrum of RX J1856.5-3754 
to NS atmosphere models of heavy-element compositions 
assuming that the NS is a fast rotator (the models are 
from ZP02). 
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Fig. 20. Spectral flux of RX 
by the best blackbody fit to 



J1856. 5-3754 as given 
the LETG data (blue 



curves, marked "hard BB"), with T{£ = 0.73 MK and 
-fi^bb = 4.4(d/120pc) km. The red curves give an ex- 
ample of a blackbody fit with T b °g s = 0.30 MK and 
g = 18.2(d/120pc) km ("soft BB") to the optical 
data (crosses). Solid curves show the absorbed spectra at 
riH = 1.0 x 10 20 cm -2 (dots are the unabsorbed models). 



a 55 ks exposure with the HRC-S/LETG instrument. A 
detailed analysis of this observation is given by Burwitz 
et al. (2001), who reported no detection of any signifi- 
cant features in the high-resolution spectrum. They also 
concluded that the combined X-ray and optical data rule 
out all the available NS atmosphere models: light-element 
models overpredict the actual optical fluxes (as first es- 
timated by Pavlov et al. 1996), whereas heavy-element 
models do not fit the LETG spectrum. In October 2001 
this source was observed for 450 ks, with the same instru- 
ment. However, even this very deep observation revealed 
no spectral features. The spectrum excellently fits a black- 
body model of T b °g = 0.73 MK and = 4.4(d/120pc) 
km^] with the hydrogen column density of Tin — 1.0 x 10 20 
cm" 2 (Fig. fji§ see also Drake et al. 2002). 

All attempts to find pulsations from J1856 have been 
unsuccessful so far (e.g., Ransom et al. 2002). However, 
since time resolution of both the ROSAT and Chandra 
data on J 1856 was not sufficient to search for periods 
shorter than » 10 ms, one could assume that the period 
of J1856 is shorter than that. Such a hypothesis might 
also explain the featureless spectrum of J1856 — spectral 
lines in radiation of a fast-rotating NS can be smeared out 



3 The distance to J1856, d ~ 120 pc, has been recently 
established by Kaplan et al. (2002a) and Walter & Lattimer 
(2002). 
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because of the Doppler effect (rotational broadening; see 
ZP02 for more details). However, atmosphere models for a 
fast-rotating NS fail to fit the detected spectrum because 
of strong deviations at energies around most prominent 
absorption edges (see Fig. and Fig. 12 in ZP02). 
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Fig. 21. Chandra LETG spectrum of RX J0720.4-3125. 
The red curve shows the best blackbody fit to the obser- 
vational data (model parameters are given in the text). 



on J1856. In this two-temperature picture, it is tempting 
to interpret the hot area as a pulsar polar cap. This seems 
to be consistent with the recent discovery of the H Q nebula 
surrounding J1856 (van Kerkwijk & Kulkarni 2001b), pre- 
sumably a bow shock in the ambient medium created by 
the supersonic motion of a NS with a relativistic (pulsar?) 
wind. However, the hot spot is unusually large and lumi- 
nous [Lg^j = 4 x 10 31 (ci/120pc) 2 erg s _1 ] for a standard 
pulsar polar cap. This could be explained by a very strong 
magnetic field localized in the spot region (i.e., essentially 
different from a centered dipole). Such a field can quench 
the pair cascade in the wind of primary particles pulled 
out from the NS surface, which explains the lack of X- 
ray/optical synchrotron radiation. On the other hand, the 
strong field can increase the heat flux from the core/crust 
to the surface. The large size of the hot spot, together 
with the GR effects (ZP02), can explain the nondetection 
of pulsations (pulsed fraction / p ^ 4% from the Chandra 
data), particularly if the hot spot is close the rotational 
pole or the direction of the rotational axis is close to the 
line of sight. If the new observation of J1856 with XMM- 
Newton (conducted in April 2002) reveals pulsations of 
this objects, it will support the two-temperature interpre- 
tation and provide constraints on the magnetic moment, 
hot spot position and inclination of the rotational axis. 



Thus, we have to admit that the blackbody spectrum 
is the best model for the X-ray spectrum of J1856. In prin- 
ciple, such a blackbody spectrum could mimic the actual 
spectrum in the relatively short band of 0.15-1.0 keV if 
the NS surface is in a condensed form (according to Lai & 
Salpeter 1997, this can happen if the magnetic field is very 
strong at the NS surface, B > 10 14 G). However, the radius 
of the emitting area, ~ 4-5 km, is too small to be the 
NS radius. Moreover, taking into consideration the opti- 
cal fluxes detected from J1856 (van Kerkwijk & Kulkarni 
2001a; Pons et al. 2002), one encounters one more prob- 
lem — the extrapolation of the blackbody model inferred 
from the X-ray data into the optical band falls well below 
the actually detected fluxes. One may assume (Pons et al. 
2002) that the X-ray radiation is emitted from a hot area 
on the NS surface, while the optical radiation is emitted 
from the rest of the surface, with a lower temperature. 
The optical spectrum is well described by the Rayleigh- 
Jeans law (van Kerkwijk & Kulkarni 2001a) and gives the 
following relation between the temperature and the ra- 
dius of the NS surface emitting the soft blackbody radi- 
ation: R^ h s ~ 8.3(T^ S /1 MK)~ 1 /2( d / 10 o pc ) km (with 
allowance for the contribution of the "hard blackbody" ) . 
On the other hand, a 3<r upper limit on the temperature 
of the soft component, TJ°£ S < 0.39 MK, found from the 
X-ray data, corresponds to a lower limit for the radius 
^bb s > 16.3 km for d — 120 pc. Figure po| shows an exam- 
ple of the two-blackbody fit to the optical and X-ray data 



. .■ -.. hard BB 




Energy (keV) 

Fig. 22. The same as in Fig. |(J for RX J0720-3125. The 
parameters for the model spectra are given in the text. The 
optical data (crosses) are from Kulkarni & van Kerkwijk 
(1998). 
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4.2. RX J0720.4-3125 

Chandra observed this object in February 2000 for 38 ks 
with the HRC-S/LETG instrument. Fitting various spec- 
tral models to these data gives results remarkably sim- 
ilar to those for J1856. Although the light-element (hy- 
drogen and helium) NS atmosphere models provide ac- 
ceptable fits (but the models are somewhat harder than 
the observed spectrum at E > 0.75 keV), the predicted 
optical fluxes (e.g., B = 23.7 and 23.3 for the hydrogen 
and helium nonmagnetic models, respectively) are much 
brighter than observed (B ~ 26.6 — see Kulkarni & van 
Kerkwijk 1998). Neither iron nor solar- mixture models 
fit the data because of numerous features in the model 
spectra. The best fit is provided by a blackbody model 
of T,°g = 0.92 MK and a radius of the emitting area of 
= 2.2 (d/lOOpc) km, absorbed with the hydrogen 
column density — 1.7 x 10 20 cm~ 2 (Fig. ^l]). How- 
ever, similar to J1856, the best blackbody model inferred 
from the X-ray data underpredicts the optical fluxes which 
can be interpreted as emitted by a "soft blackbody" with 
R^ a ~ 4.3 (T b °g s /1 MK)- 1 / 2 (d/100pc) km. To obtain an 
acceptable fit of the two-blackbody model to the X-ray 
data, one should keep Tgg s ^ 0.50 MK, that results in 
^bb s ~ 6.1 (rf/100pc) km. Figure shows an example 
of the two-blackbody fit to the optical and X-ray data, 
with T b °g s = 0.45 MK and R^ h s = 6.4(d/100pc) km for 
the soft component. If future optical observations provide 
the distance to J0720. 4-3125, it will be possible to eval- 
uate the NS radius and, based on the X-ray light curves, 
estimate the mass-to-radius ratio. 

5. Concluding remarks 

Because of the limited space, we described the observa- 
tional results on only a handful of observed thermally 
emitting isolated NSs. We only briefly mentioned AXPs 
and SGRs, whose X-ray radiation almost certainly con- 
tains a thermal component, with temperatures of 5-10 
MK and equivalent radii of 1-5 km. The current status of 
AXPs, including their comparison with SGRs, is reviewed 
by S. Mereghetti et al. (this volume). From the observa- 
tional perspective, the nature of these enigmatic objects 
can hardly be understood without deep multiwavelength 
observations, particularly in the IR and hard X-ray bands, 
which would help elucidate the origin of the nonthermal 
component and separate the thermal component more pre- 
cisely. We did not touch "elderly" nearby pulsars (e.g., 
B1929-10, B0950+08), mostly because they have not been 
observed with Chandra. If future deep observations show 
that their X-ray radiation is indeed thermal (emitted from 
hot, small polar caps), as suggested by the ROSAT and 
ASCA observations, the analysis of the spectra and light 
curves will help establish the strength and geometry of 
their magnetic fields and the NS mass-to-radius ratio. We 
did not mention important non-detections of thermal ra- 
diation from very young pulsars, the Crab (Tennant et al. 



2001; M.C. Weisskopf, this volume) and the recently dis- 
covered PSR J0205+6449 in the SNR 3C 58 (Slane et al. 
2002), which put upper limits on the NS surface tempera- 
ture. Finally, very interesting results on thermal radiation 
from transiently accreting NSs in quiescence have been 
reported recently (Rutledge et al. 2001, and references 
therein). Although these old NSs are not truly isolated, 
their radiation in quiescence resembles that of young iso- 
lated NSs because of additional heating due to pyenonu- 
clear reactions in the compressed accreted material. 

The Chandra observations described above have shown 
a few important things. First of all, they left no doubts 
that thermal emission indeed dominates in the soft X-ray 
radiation of quite a few isolated NSs, and the analysis 
of this radiation allows one to measure the NS temper- 
atures^] and, in some cases, the radii or mass-to-radius 
ratios. Particularly important is the conclusion that the 
simple picture of a NS with a centered dipole magnetic 
field and uniform surface temperature is, most likely, an 
oversimplification. We see that whenever the data allow 
a detailed analysis, the temperature is not uniform, in 
both active pulsars and radio-quiet NSs. Moreover, the 
radius of X-ray emitting area is, as a rule, considerably 
smaller than the "canonical" NS radius. The example of 
J1210-5226 (Sec. 3.1) shows that the characteristic age 
(t c = P/2P) of a pulsar can differ from its true age by a 
large factor, and the conventional "pulsar magnetic field" 
can be quite different from the actual magnetic field at 
the NS surface. Therefore, any inferences on the proper- 
ties of NSs and the superdense matter obtained without 
taking this into account should be considered with cau- 
tion. For example, adopting the X-ray emitting radius to 
be the NS radius, a number of authors hurried to conclude 
that NSs are strange quark stars. Another example is the 
numerous comparisons of the NS cooling curves with ob- 
servations of thermal X-rays from pulsars postulating that 
the pulsar's characteristic age is the true age and the em- 
pirical temperature of one of the "thermal components" 
is a single effective temperature. If we want to understand 
the real objects, we should abandon the naive paradigm 
and answer the basic questions raised by the observations: 

• What is the true explanation of the small, hot ar- 
eas in many, if not all, NSs, such as J2323+5848 
(CCO of Cas A), RX J1856.5-3754, etc? If these small 
heated areas are associated not with (unseen) pul- 
sar activity, but with superstrong local magnetic fields, 
how such fields affect the thermal evolution of NSs? 

• What is the relation between the various kinds of NSs 
- CCOs, AXPs, SGRs, active pulsars, dim NSs? For in- 
stance, why do CCOs, AXPs and SGRs show similar prop- 
erties of thermal radiation, being so different in the other 
observational manifestations? Can it be that the mem- 



4 See Fig. 5 in Mereghetti et al. (this volume) which demon- 
strates current estimates for the temperatures and luminosities 
of 22 thermally emitting NSs. 
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bers of different subclasses are intrinsically similar objects 
at different stages of their evolution (e.g., AXPs are de- 
scendants of CCOs) or viewed from different directions 
(e.g., dim NSs are usual middle-aged pulsars whose pulsar 
radiation is not seen because of unfavorable geometry)? 

• What is the actual temperature distribution over 
the NS surface? Does the two-component model for 
thermal radiation we have used (thermal soft + 
thermal hard) describe the distribution adequately? 

• What is the actual mechanism(s) of the surface emis- 
sion? Are there any gaseous atmospheres at all? Why does 
the blackbody model describes the observed spectra so 
well in many cases? 

These are the questions to be addressed (and, hopefully, 
answered) by future observations of NSs. 
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